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ABSTRACT 

We have obtained ~ 100 pc resolution CO and ~ 60 pc resolution Ha observations of the dwarf 
spiral galaxy NGC 4605. We use them to derive a high resolution rotation curve and study the 
central density profile of NGC 4605's dark matter halo. We find that these observations do not 
agree with the predictions of most high resolution Cold Dark Matter calculations. We investigate 
two extreme cases: 1) NGC 4605 has a maximal exponential disk, that we model using if-band 
observations and remove to study the structure of its dark matter halo, and 2) NGC 4605 is dark 
matter dominated and its disk is dynamically negligible. Because the mass-to-light ratio of the 
maximal disk is already very low we favor the first solution, which indicates the halo has one 
component with a density profile p cx r~ - 65 out to R ~ 2.8 kpc. In the second case, the rotation 
curve requires the presence of two components: a small (~ 600 pc) core surrounded by a much 
steeper p oc r -11 halo. Removal of intermediate (submaximal) disks does not ameliorate the 
discrepancy between the predictions and the observations. 



1. Introduction 

The Cold Dark Matter (CDM) paradigm and 
its variations (e.g., ACDM) have enjoyed remark- 
able success at explaining the large scale struc- 
ture of the Universe. At the scale of individual 
galaxies, however, CDM faces several serious chal- 
lenges (e.g., Sellwood & Kosowski 2000). One of 
the key predictions of CDM and ACDM simula- 
tions is that all dark matter halos share a uni- 
versal density profile (Navarro, Frenk, & White 
1996; Huss, Jain, & Steinmetz 1999). The origi- 
nal calculations predicted a central density profile 
p ~ r" 1 , where r is the galactocentric distance. 
More recent, higher resolution simulations have 
found even steeper p ~ r~ 15 central cusps (e.g., 
Moore et al. 1999; Ghigna et al. 2000). 

These calculations, however, collide with a 
number of observational studies of dark-matter 
density profiles. Most of the observations find ei- 
ther central constant density cores (e.g., Flores 
& Primack 1994; McGaugh & de Blok 1998; de 
Blok et al. 2001; Bordello & Salucci 2001), or pro- 
files where the density increases more slowly than 
r _1 toward the center (e.g., Kravtsov et al. 1998). 



Other authors conclude that the observational evi- 
dence available to date is inconclusive, largely due 
to the lack of high quality high resolution rota- 
tion curves (van den Bosch et al. 2000). Indeed, 
using H I interferometers to obtain rotation curve 
data with the angular resolution and sensitivity 
required in order to resolve the central cusps is a 
daunting task for all but the nearest galaxies. The 
usual alternative, deriving rotation curves from 
longslit Ha spectra, provides spatial information 
only along one dimension. This makes it impossi- 
ble to separate circular and non circular motions 
and to check for position angle errors. Other po- 
tential problems with rotation curves based solely 
on Ha are the systematic effects introduced by 
extinction and the local expansion flow of H II 
regions, which are especially important close to 
galaxy centers, precisely where the shape of the 
rotation curve is most important. A recent study 
by Blais-Ouellette, Amram, & Carignan (2001), 
however, shows good agreement between the mea- 
sured H I and Ha kinematics in the inner parts 
of their galaxies, indicating that at least in some 
galaxies extinction is not a substantial concern. 
Despite the recent body of work on the subject, 
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it remains clear that we need more high quality, 
high resolution rotation curve measurements to 
advance towards solving the galactic scale dark 
matter puzzle. 

Since both H I and Ha rotation curve studies 
are problematic, we have opted for a third, com- 
plementary technique: we use CO interferometry 
to measure the inner rotation curves. Millime- 
ter wave interferometry avoids the problems inher- 
ent to the first two techniques. First, it is much 
more sensitive to surface brightness than H I in- 
terferometry, allowing us to routinely obtain sen- 
sitive data with 5" angular resolution or better 
while the typical resolution of most H I studies 
is ~ 30". Second, not only is CO interferometry 
unaffected by extinction, but it works best near 
galaxy centers where most molecular clouds are 
typically found. Third, like H I interferometry but 
unlike long slit Ha spectroscopy, it provides com- 
plete spectral imaging information allowing us to 
obtain reliable rotation curves. 

The relative sensitivity advantage of CO in- 
terferometry stems from the fact that the 12 CO 
(J = 1 — > 0) rotational transition at 2.6 mm has 
a wavelength ~ 80 times smaller than the 21 cm 
spin flip H I transition. This allows the CO an- 
tennas to be placed <~ 80 times closer together 
than the H I antennas for the same angular res- 
olution. Because the sensitivity of aperture syn- 
thesis to surface brightness is proportional to the 
fraction of the filled aperture, compactness confers 
a considerable sensitivity advantage to millimeter 
wave interferometers. For example, the ten ele- 
ment Bcrkcley-Illinois-Maryland Association mil- 
limeter array (BIMA; Welch et al. 1996) has a 250 
times larger filled aperture fraction than the VLA. 

The difficulty in using CO interferometry for 
dark matter studies, however, is finding suitable 
target galaxies. Late type dwarf galaxies are ideal 
for such studies, because their density profiles are 
dominated by dark matter in their inner regions 
(Persic & Salucci 1990; Salucci, Ashman, & Persic 
1991). Although little is known about the molec- 
ular gas content of most dwarf galaxies, with a 
few notable exceptions they are not strong molec- 
ular sources. We have obtained CO interferometer 
maps of several nearby dwarf galaxies with single- 
dish detections, a subset of which we selected for 
rotation curve studies. To further expand the set 
of suitable sources we are carrying out a system- 



atic survey of nearby dwarf galaxies in CO. Thus, 
we expect that the results presented here are only 
a first glimpse of the final study. We present our 
observations in §2, perform a rotation curve study 
in §3, and summarize our results in §4. 

2. Observations and Data Processing 

The source we discuss in this work, NGC 4605, 
is an SBc dwarf galaxy. Independent photometric 
estimates place this galaxy ~ 4.6 Mpc (Karachcnt- 
sev & Tikhonov 1994) and 4.26 ± 0.64 Mpc away 
(M. Pierce, priv. comm.), yielding a physical scale 
of ~ 20 pc/". The physical resolution of our CO 
observations (and of our rotation curve) is there- 
fore ~ 100 pc. The extinction corrected abso- 
lute B magnitude of NGC 4605 is M B = -18.10 
and its size is R25 ~ 3' (~ 3.6 kpc), very similar 
to the magnitude and size of the Large Magel- 
lanic Cloud (LMC). The line-of-sight inclination 
of NGC 4605 is i w 69 deg, and H I observations 
find an inclination corrected linewidth at 20% in- 
tensity W 20 « 198 km s -1 (Bottinelli et al. 1985). 
This indicates a maximum rotational velocity of 
~ 100 km s -1 (again, similar to that of the LMC) 
that along with its size and luminosity definitively 
identifies this galaxy as a dwarf. Carbon monoxide 
millimeter emission was detected in NGC 4605 by 
Young et al. (1995) using the FCRAO 14 m radio 
telescope. Our first interferometer maps revealed 
CO emission throughout the disk of NGC 4605, 
which became the obvious first choice for a rota- 
tion curve study. 

We observed the 12 CO ( J = 1 — > 0) rotational 
transition in NGC 4605 using three configurations 
of the BIMA array (B, C, and D) between June 
2000 and March 2001. We mapped one BIMA pri- 
mary beam, subtending a field with a half power 
diameter of ~ 100", sufficient to encompass the 
CO emission. The spectrometer was configured 
with 2 km s -1 wide channels and a bandpass 
of 260 km s -1 . The individual tracks were cali- 
brated, combined, imaged, and deconvolved using 
the CLEAN algorithm within the MIRIAD astro- 
nomical package. The resulting naturally weighted 
map has a beam size of 4.8" x 5.4" (~ 99 x 111 
pc) with P. A. « 6 deg (Fig. 1). The individual 
planes of the datacube have an RMS of ~ 24 mJy 
beam -1 (~ 85 mK) at a velocity resolution of 3 
km s . 
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The datacube was used to produce a first mo- 
ment map of the CO emission (Fig. 2). The 
map was rotated to align the optical major axis 
of the galaxy with one of the coordinate axes 
(P. A. — 125 deg), and deprojected using the opti- 
cal center and inclination angle. The velocity data 
were corrected by the sine of inclination angle, and 
sorted in concentric rings of width 5". To each ring 
we fit a model v t, s = v t + v c cos + v r sin 6, rep- 
resenting the projected effects of translation (v t ), 
circular (v c ), and radial (v r ) velocities. The er- 
rors in the velocity map, multiplied by the square 
root of the beam oversampling factor, are used to 
weight each point in the fit (Fig. 3). The er- 
ror bars for the fits (which range from 1.4 to 3.6 
km s _1 ) are computed using the diagonal of the 
covariance matrix. We find a small v r component, 
increasing outwards to ~ 25" and then leveling 
off (Fig. 4). This component is never larger than 
half v c , and adding it in quadrature to v c does 
not significantly change the density profile expo- 
nents discussed in this paper. The presence of this 
component may indicate that the kinematics in- 
side ~ 25" have a position angle slightly different 
from the optical P. A. (closer to P. A. ~ 95 deg). 

To extend the CO rotation curve beyond 50" 
(~ 1 kpc), we observed NGC 4605 with the 0.6 m 
Coude Auxiliary Telescope at Lick Observatory on 
the night of 23 April 2001. We used the Hamil- 
ton Echelle Spectrometer with a 75A wide filter 
centered on Ha (A = 6562. 8A) to block adjacent 
spectral orders. The slit was 640 /mi (6") wide 
and ~ 8' long, oriented along NGC 4605's major 
axis at P. A. = 125 deg using the facility's focal 
plane image rotator. Our velocity resolution was 
~ 6 km s _1 , with a spatial resolution of ~ 3". 
This spectral resolution is much higher than that 
of most published rotation curve studies, and it 
allows us to very accurately measure the kinemat- 
ics of the inner disk of NGC 4605. We integrated 
for 4.5 hours divided in 30 minute exposures using 
manual guiding. 

We reduced the Ha data using IRAF. We re- 
moved a bias frame from the individual exposures, 
and combined them using the imcombine rou- 
tine with its built in cosmic ray rejection algo- 
rithm. Because no continuum was detected in 
the individual images, we registered them using 
offsets obtained by cross correlating the spectra. 
For wavelength calibration we used exposures of 



a ThAr reference lamp, and observations of a ra- 
dial velocity standard star. To extract the rota- 
tion velocities we fit a gaussian to each row of the 
spectrum. We folded the results into a rotation 
curve by defining the center of the galaxy to be 
the point at which the two sides of the rotation 
curve lined up best. The typical errors in the in- 
dividual points of the folded rotation curve were 
less than 2 km s _1 . Comparison of the CO and Ha 
rotation data reveals that both tracers are in good 
agreement over the common region (Fig. 4). Thus 
Ha appears to be a good tracer of the kinematics 
of the inner disk in NGC 4605. 

We developed a simulation to estimate how our 
finite angular resolution affects both the CO and 
Ha rotation curves (beam smearing). Because 
calculating the effect of the beam smearing re- 
quires assuming an underlying velocity field, we 
performed these computations in an iterative, self 
consistent manner. The net effect is very small, 
amounting to a velocity loss of ~ 2 km s _1 for the 
innermost points and rapidly decreasing outwards. 
Although small, this systematic velocity loss is not 
entirely negligible in the inner 30". We removed 
the effect of beam smearing from the data before 
performing the fits. 

We used 2MASS K s data to derive the distri- 
bution of mass in the disk of NGC 4605. The 
advantage of using K -band observations to model 
the disk mass over more typical B, V, or / data 
cannot be overemphasized: at if-band the extinc- 
tion is a factor of ~ 10 down from V, and a fac- 
tor of ~ 5 down from I (Rieke & Lebofsky 1985). 
We retrieved the archival 2MASS atlas image for 
NGC 4605 and extracted the surface brightness 
profile using the photometric calibration provided 
in the FITS header. We then fit an exponential 
disk to the surface brightness profile using data 
with galactocentric radii r G [10", 60"], finding 
[iK ~ 16.35 + 0.029 r where is the apparent 
-fT-band magnitude. NGC 4605 does not have a 
significant bulge, and the surface brightness distri- 
bution is well fit by an exponential law over this 
range. 

3. Results and Discussion 

The universal dark matter halo density law pro- 
posed by Navarro, Frenk, & White (1995) follows 
the form 
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P( r ) = S c n , 

Peru (r/r s ){l + r/r s f ' [> 

where p crlt = 3H 2 /8ttG ~ 1 x 10~ 29 g cm~ 3 is the 
critical density, 5 C is the halo overdensity, and r s is 
the characteristic radius (simulations suggest r s of 
several kpc). At galactocentric distances r <C r s 
the halo has an approximate density profile p w 
Per a 5 el's r . In general, a density profile p oc 
r~" implies 

Vc oc r< 2 -")/ 2 . (2) 

Most dark matter density profile studies assume 
that the halo is the dominant mass component, 
and do not attempt to model and remove the con- 
tribution from the galactic disk to the rotation, 
which is taken to be negligible. We call this ap- 
proach the "minimal disk" hypothesis, since it im- 
plies that the mass-to-light (M/L) ratio of the 
material in the disk is M/L ~ 0. The opposite 
limit is given by the "maximal disk"; the maxi- 
mum exponential disk that can be subtracted from 
the measured rotation curve without creating neg- 
ative residual velocities. Fitting a maximal disk 
to the rotation curve yields a maximum mass-to- 
light ratio for the material in the disk (predomi- 
nantly stars and gas). 

To find the upper limit to the M/L ratio in 
NGC 4605 we fit a thin maximal exponential disk 
using our rotation curve and the 2MASS if-band 
data. The fit is constrained by the innermost 
points of the rotation measurements. This fit 
yields M/L K = 0.25 ± 0.03 M Q /L QK (we as- 
sume Mqk = 3.34), implying a disk mass surface 
density of stars plus gas p = 0.053 exp(— r/37'.'l) 
g cm -2 . The extinction corrected color, B — R = 
0.67 (M. Pierce, priv. comm.), predicts a simi- 
lar mass-to-light ratio M/L K ~ 0.34 (Bell & dc 
Jong 2001). This value is very close to the lower 
limit of the reasonable range of M/Lk, possibly 
corresponding to a young starburst population ac- 
cording to the Starburst99 models (Leitherer et al. 
1999). This idea is further supported by the obser- 
vation that, in dwarf galaxies, bright CO emission 
appears to be accompanied by the intense star for- 
mation activity signaling a young starburst. The 
fact that the upper limit of M/L is so low suggests 
that the disk of NGC 4605 is quite close to max- 
imal, as submaximal disks will imply even lower 



values of the mass-to-light ratio. Using the sin- 
gle dish CO and H I data and allowing for Mag- 
ellanic metallicities we estimate the gas is ~ 20% 
of the total mass of the maximal disk. The effect 
of including the gas would be to lower the stellar 
mass-to-light ratio to M/Lk ~ 0.20 Mq/Lqk- 
Because this correction is somewhat uncertain we 
ignore it in the following discussion, but note that 
its only effect is to make the M/L ratio even lower. 

The rotation due to an exponential disk is 
straightforward to compute (Freeman 1970), and 
to remove it from the measured rotation we use 
v c = v Lio + v diski where a11 velocities are cir- 
cular velocities (Binney & Tremaine 1987). Fig- 
ure 5 shows the result of removing the maximal 
disk contribution from the overall rotation. To 
the residual rotation, due to the halo, we fit a 
Vhalo K rl3 power law in the region r G [40", 140"] 
and find the corresponding a using Eq. 2. The 
best fit power law is v c ~ 3. Or 067 km s _1 , im- 
plying p ~ 4.9 x 10- 24 i?~ - 65 g cur 3 (where R 
is the galactocentric radius in kpc). This value 
of a is significantly lower than the calculations by 
Navarro et al. (1996) or Moore et al. (1999) in- 
dicate. Note, however, that it agrees well with 
the critical solution recently found by Taylor & 
Navarro (2001) with a ~ 0.75. Thus, removing the 
maximal disk from the rotation curve of NGC 4605 
results in a shallow halo, with a density exponent 
a = 0.65 ± 0.02. The error bars including the un- 
certainty in the M/L ratio are ±0.07. 

It is interesting to investigate whether the halo 
rotation of NGC 4605 shows any hint of a con- 
stant density core, such as those proposed for other 
galaxies (e.g., Borriello & Salucci 2001). The re- 
gion inward of 30" in Fig. 5b does not show a 
significant departure from the outer halo fit, es- 
pecially after correcting the innermost points for 
the effects of beam smearing. While most veloc- 
ity residuals inside 30" are negative, the overall 
deficit is less than 2 km s _1 . Because the maxi- 
mal disk is dynamically dominant inside r = 40" 
any departures from a perfect exponential disk, or 
from a constant mass-to-light ratio, can account 
for it. Among other problems that could intro- 
duce similar systematic errors near the centers of 
galaxies are: 1) a displacement of the Ha slit by a 
few arcseconds perpendicular to the major axis of 
the galaxy, 2) a disagreement of a few arcseconds 
between the optical and dynamical centers, or 3) 
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the presence of a weak galactic bar with the conse- 
quent non circular orbits. For these reasons, it is 
difficult to completely exclude a constant density 
core using these data: in the presence of a dynam- 
ically dominant disk the rotational signature of 
such core is a deficit of a few km s _1 throughout 
this region. Thus, although we find no evidence 
for a constant density core, we cannot discard it 
completely. 

What happens if despite the evidence for a max- 
imal disk in NGC 4605 we assume the minimal 
disk hypothesis? It is easy to see that the observed 
rotation cannot be fit well by a single power law 
(Fig. 4). This suggests that there are at least 
two distinct mass components that are contribut- 
ing to the potential. Indeed, a possible decom- 
position of the rotation curve of NGC 4605 indi- 
cates a slowly varying density profile inside 30" 
(i.e., a core with v c ~ 3.0 r 8 km s _1 , implying 
/) ~ 1.5 x 10~ 23 R~ 0A g cm -3 ), and a rapidly de- 
creasing density outwards of 40" (v c ~ 10.4 r 45 
km s _1 , yielding p~8x 10~ 24 R^ 11 g cm -3 ). In 
the absence of a disk the observed rotation requires 
a two component halo: a core with a = 0.4±0.2 in- 
side r ~ 30" and an outer halo with a = 1.1±0.02. 

Both of these limits are in disagreement with 
most results of collisionless cold dark matter sim- 
ulations. The removal of a submaximal disk does 
not eliminate the discrepancy Such a solution has 
a combination of the problems of each limit, re- 
quiring both a shallow outer halo and an inner 
core. It should be noted that the goodness of fit 
of the outer halo by a power law (or any other 
mathematical model used to describe it) cannot 
be reliably employed to obtain a "best fit" mass- 
to-light ratio and select a solution. Indeed, we 
find that for M/L in the range M/L e [0, 0.25] 
Mq/Lqk, the RMS of the residuals of the power 
law fit to the rotation curve only changes from 1.5 
(Fig. 4b) to 1.9 km s _1 (Fig. 5b), showing that the 
RMS is very insensitive to how much of the stellar 
disk is removed. We have argued that, because we 
find a low mass-to-light ratio, NGC 4605's disk is 
close to maximal. This solution has the virtue of 
requiring only one halo component. Regardless of 
the precise M/L ratio of the disk of NGC 4605, 
the inescapable conclusion is that the shape of the 
halo of NGC 4605 poses a challenge to the CDM 
simulations. 



4. Summary and Conclusions 

We have measured the kinematics of the dwarf 
galaxy NGC 4605 with high spectral (~ 2 km s _1 
for the CO, ~ 6 km s _1 for the Ha) and spa- 
tial (~ 100 pc) resolution. Our study of the com- 
bined CO and Ha rotation curve leads to a ro- 
bust conclusion: the standard collisionless CDM 
simulations do not accurately predict its shape. 
We have discussed two extreme cases: 1) the mea- 
sured rotation is due to a combination of the halo 
and an exponential maximal disk, and 2) the ro- 
tation is due to solely the halo (i.e., the minimal 
disk hypothesis). Because the if-band mass-to- 
light ratio we find for the first case is very low 
(M/L K ps 0.25 Mq/Lq K ), we believe that the 
disk of NGC 4605 is close to maximal and that 
the first solution is the correct one. This solution 
requires only one halo component within R ~ 2.8 
kpc. The resulting p ~ r~ a halo density profile 
has a = 0.65±0.07. This exponent is substantially 
lower than the predictions of most high resolution 
CDM calculations, although it agrees much bet- 
ter with the critical solution found by Taylor & 
Navarro (2001). 

If we disregard the evidence favoring a maxi- 
mal disk and adopt the minimal disk case, we find 
that we need a soft p ~ r~ 0,4 core to explain the 
rotation of the innermost 600 pc. Note that while 
the precise value of a for this core is rather uncer- 
tain, the fact that a break in a is required is not. 
Furthermore, the density profile of the halo in the 
0.8 to 2.8 kpc region has a = 1.1 ± 0.02. This is 
still noticeably different from the results of some of 
the recent high resolution simulations which find 
a ps 1.5 (Moore et al. 1999; Ghigna et al. 2000), 
although it agrees with other calculations that find 
a ps 1.1 (Jing & Suto 2000). None of these calcu- 
lations, however, predict the core. The removal of 
a submaximal disk has a combination of the prob- 
lems of these two extreme cases. 

Notwithstanding which solution is the best, 
both of these extremes (and all of the intermediate 
cases) are discrepant with most of the available 
predictions of standard CDM. We have to note 
that the strength of these conclusions stems from 
the fact that both the CO and the Ha data, in- 
dependently, reveal the same kinematics for the 
inner 50" of NGC 4605. A similar analysis, rely- 
ing only on one or the other rotation curves, would 
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be impaired by the possible systematics of the in- 
dividual datasets. 
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Fig. 1 — 2MASS if-band image of NGC 4605 
with integrated intensity CO contours overlaid. 
The RMS of the integrated intensity data is ~ 130 
mJy beam -1 km s" 1 (~ 450 mK km s" 1 ). The 
contours start at 3a and increase in steps of 4cr to 
3 Jy beam -1 km s" 1 . The offsets are with respect 
to NGC 4605's center, at a J2 ooo = 12 h 40 m 00!0, 
<Jj20oo = 61°36'31". The synthesized beam is 
shown in the upper left corner. 
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Fig. 3. — Individual fits to the CO rings, velocity 
plotted against angle. The number in the upper 
left corner indicates the ring number, increasing 
outwards. The vertical lines are the individual 
data points with their errors, including the square 
root of oversampling factor. The gray thick line is 
the overall fit to w h s , while the black thin line is 
the v c component, as described in the main text. 



Fig. 2. — Velocity field obtained from the CO dat- 
acube. The typical la error in the velocities is 
^2 — 3 km s _1 . The ellipses correspond to the 
projection of the 5" wide rings used to determine 
the rotation curve. 
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Fig. 4. — Combined CO and Ha rotation curve of 
NGC 4605 with no disk removed, in linear-linear 
scale. Plot (a) shows the individual rotation mea- 
surements for CO (filled squares) and Ha (open 
squares). It also shows the v r component obtained 
from the CO (filled triangles). The CO is sampled 
every 5", the Ha every 3", and both have been 
corrected for the effects of beam smearing. The 
solid line is the power law fit to the outer halo Ha 
data in the r G [40", 140"] region, p oc r~ 1A . The 
dashed line is the power law fit to the CO rotation, 
yielding p oc r~ - 6 . The gray thick line is the max- 
imal disk obtained from the if -band observations. 
Plot (b) shows the residuals in the data and the 
fits, after removing the outer halo fit. The dotted 
line is the effect of beam smearing, as it was re- 
moved from the Ha observations. The correction 
applied to the CO data is very similar. The error 
bars are internal, statistical errors obtained from 
the fits. 
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Fig. 5. — Combined CO and Ha rotation curve of 
NGC 4605 after removing the maximal disk. The 
key to symbols and lines is the same as in Fig. 4. 
The few data points missing had v c < Vdisk , yield- 
ing imaginary v^aio- Plot (a) shows the new outer 
halo fit for r G [40", 140"], p oc r -°- 65 , and the 
fact that the maximal disk is dynamically domi- 
nant inside r ~ 50". Plot (b) shows the residuals 
after removing the outer halo fit. The gray re- 
gions indicate la, 2a, and 3a departures from the 
fit. There is no evidence for a break in the power 
law (i.e., a core) after the disk is removed. 
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